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MINERALOGY 


Discovery of davemaoite, CaSi03-perovskite, 
as a mineral from the lower mantle 


Oliver Tschauner"*, Shichun Huang’, Shuying Yang”, Munir Humayun’, Wenjun Liu’, 
Stephanie N. Gilbert Corder*, Hans A. Bechtel*, Jon Tischler*, George R. Rossman? 


Calcium silicate perovskite, CaSi03, is arguably the most geochemically important phase in the lower 
mantle, because it concentrates elements that are incompatible in the upper mantle, including the 
heat-generating elements thorium and uranium, which have half-lives longer than the geologic history of 
Earth. We report CaSi03-perovskite as an approved mineral (IMA2020-012a) with the name davemaoite. 
The natural specimen of davemaoite proves the existence of compositional heterogeneity within the 

lower mantle. Our observations indicate that davemaoite also hosts potassium in addition to uranium 
and thorium in its structure. Hence, the regional and global abundances of davemaoite influence the 
heat budget of the deep mantle, where the mineral is thermodynamically stable. 


alcium silicate, CaSiO3, occurs in a variety 

of natural and synthetic polymorphs. The 

low-pressure polymorph wollastonite is 

a common metamorphic mineral. Breyite 

(J), an intermediate-pressure polymorph 
(2), has been found as inclusions in diamond. 
At the pressures and depth range of Earth’s 
transition zone (420 to 660 km) and lower man- 
tle (LM; 660 to ~2700 km), CaSiO; assumes a 
perovskite structure. Perovskite-type CaSiOs, 
first synthesized by Liu and Ringwood (3), is a 
liquidus phase for basaltic and peridotic bulk 
rock compositions at LM pressures and tem- 
peratures and has been experimentally shown 
to host many elements that are incompatible 
in upper-mantle minerals (4-7). These include 
rare-earth elements (REEs), large ion lithophile 
elements (LILEs; K, Sr, and Ba), Ti, U, and Th. 
In other words, these elements are compatible 
rather than incompatible in a LM mineral as- 
semblage that contains a few vol % of CaSiO;- 
perovskite. The composition and abundance 
of this phase in the LM are therefore key in 
constraining the budget and distribution of 
REEs and LILEs and the elements with abun- 
dant radioactive isotopes (K, U, and Th) that 
make an important contribution to the heat 
of Earth’s mantle (8). Through these param- 
eters, CaSiO; perovskite provides essential con- 
straints on the fate of recycled crust in deep 
Earth, thermochemical anomalies, and the 
existence of a magma ocean at the base of 
Earth’s mantle. The synthetic perovskite phase 
of pure CaSiO; has been found to assume either 
cubic or tetragonal symmetry (9) and belongs 
to the tausonite (SrTiO3)-type perovskites that 
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adhere to fundamentally different structural 
distortion mechanisms (JO) and crystal-chemical 
constraints than the GdFeO;-type perovskites 
such as bridgmanite [MgSiO;-perovskite (77)] and 
the CaTiO; mineral, actually named perovskite. 

The difficulty of finding CaSiO;-perovskite 
in nature stems from its stability at pressures 
only above 20 GPa (3, 4) along with a low ki- 
netic barrier for back-conversion into low- 
pressure polymorphs (2). This barrier is lower 
than for bridgmanite, which has been found as 
arare occurrence in highly shocked meteorites 
despite its stability only above 23 GPa (J/, 12). 
Nestola et al. (13) reported the presence of 
CaSiO3-perovskite as an inclusion in a diamond 
from the Cullinan mine, South Africa. The re- 
ported phase deviates from synthetic CaSiOs- 
perovskite in several ways: (i) its volume at 
ambient conditions is >20% larger (9); (ii) it 
sustains the beam of an electron microscope, 
whereas any synthetic CaSiO;-perovskite vitri- 
fies rapidly at ambient conditions; (iii) its cell 
axis ratios and Raman spectrum are nearly 
equal to those of CaTiO3; and (iv) its space 
group indicates a structural distortion mech- 
anism different from that of synthetic CaSiO;- 
perovskite but much closer to that of CaTiO; 
(7, 8). Nestola et al. (13) proposed that this 
distinctive phase of CaSiO; is the result of 
partial decomposition of a Ti-bearing CaSiOs;- 
perovskite. The coexistence of CaTiO; + CaSiO; 
polymorphs in diamond inclusions may also 
point to retrograde transformation of stoichi- 
ometric Ca-Si-Ti-perovskites (1) that form in 
the deep upper mantle (5 to 10 GPa). 

The findings by Nestola et al. (13) are no- 
table by themselves but differ from the ex- 
pected high-pressure CaSiO-perovskite and 
have not resulted in the approval of CaSiO;- 
perovskite as a mineral. We report the dis- 
covery of CaSiO3-perovskite as a mineral ap- 
proved by the Commission of New Minerals, 
Nomenclature and Classification (CNMNC) 
of the International Mineralogical Association 
(IMA). The new mineral [IMA2020-012a (14)] 
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is named “davemaoite” in honor of Dave (Ho- 
kwang) Mao for his eminent contributions 
to the field of deep-mantle geophysics and 
petrology. The type material—inclusions in a 
diamond from Orapa, Botswana—is deposited 
in the Natural History Museum Los Angeles 
(catalog number NHMLA 74541, formerly 
GRRI1507 of the Caltech mineral collection) 
(15). Davemaoite coexists with orthorhombic 
carbonaceous a-iron and wiistite (Feo sMgo2)O 
at 8 to 9 GPa remnant pressure (Fig. 1). Sep- 
arate inclusions of ilmenite, iron, and ice-VI 
in the same diamond (/6) have remnant pres- 
sures of 7 GPa and 8 to 9 GPa, respectively. The 
x-ray diffraction (XRD) pattern of davemaoite 
is that of a cubic perovskite (Fig. 1) (75), with, 
at most, contributions of <5 vol % of material 
with +2.5% smaller or larger volume, whereas 
an overall distortion of the lattice can be ex- 
cluded on the basis of the reflection inten- 
sities (Fig. 1). Cubic ABO; perovskites have no 
internal structural degrees of freedom and 
comprise only one chemical formula unit. Thus, 
the identification of this phase is unambiguous 
even in diffraction patterns with contributions 
from more than one phase (Fig. 1). 

Davemaoite was identified through the XRD 
pattern of cubic perovskite at a location in the 
hosting diamond with a CaKa x-ray fluores- 
cence (XRF) signal far above background. Both, 
XRF and XRD data were obtained at beamline 
34-ID-E at the Advanced Photon Source (J5). 
We superimposed the CaKo. XRF map (Fig. 2) 
on a visible light image of the holotype mate- 
rial at the beamline where we subsequently 
collected the XRF and XRD data. We also made 
corresponding maps of Fe and Ti (fig. S1). Areas 
with XRF signal at noise level in Fig. 2 show 
no x-ray diffraction besides that of diamond. 
Right after acquisition of the XRF map, we 
examined by XRD the inclusions that were 
found by XRF. XRD and XRF were collected 
on the inclusions when they were fully en- 
trapped in the doubly polished platelet of the 
hosting diamond. We focused the x-ray beam 
to an area of 0.5 um by 0.5 um in order to iden- 
tify inclusions with high spatial resolution. The 
inclusions of 4 um by 6 um and 4 um by 16 um 
areas within the red circle of the Ca XRF map 
corresponded to the XRD patterns of the cubic 
perovskite (Fig. 1). We added frames with 
perovskite patterns to obtain better signal 
and powder statistics. 

We confirmed the identification of davemaoite 
by infrared spectroscopy. Cubic ABO; perov- 
skites have no Raman-active modes and three 
infrared (IR)-active modes (15). We observed 
two of the three IR-active modes (Fig. 1, inset), 
while the third one was below the diffraction 
limit for objects as small as these inclusions. 
As expected, we observed no Raman peaks 
(fig. S2). We calculated mode energies by fit- 
ting force constants to match ab initio cal- 
culated zone-center phonon energies of the 
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tetragonal structure (17) and mapped the 
tetragonal Brillouin zone onto that of the 
cubic structure. Intensities are based on the cal- 
culated phonon density of state. 

Subsequently, we used laser-ablation induc- 
tively coupled mass spectrometry (LA-ICP-MS) 
to excavate and analyze the chemical composi- 
tions of two of the inclusions with a 100-1m- 
diameter laser beam. We indicate the ablation 
area with a red circle (Fig. 2). We monitored 
all mass peaks under medium mass resolution 
(m/Am = 4000), which allowed us to resolve 
numerous carbon-related molecular interfer- 
ences on low-mass isotopes. We hit two in- 
clusions of davemaoite at 5 to 8 um and 80 to 
100 um depth below the polished surface (Fig. 2, 
inset). The time-resolved “Ca signal of the LA- 
ICP-MS measurement (Fig. 2, inset) shows 
signal clearly above background level. We also 
collected equivalent profiles for °°Fe, ?°K, and 
52Cr (fig. SI). In all cases, the signals rose above 
background at the same depth, consistent with 
their origin in the same two inclusions. We ob- 
tained an average davemaoite composition of 
(Cao.43a)Ko.200yNao.o6F€o.1ayAlo.osM80.06CT0.04(2)) 
(Siz.o@)Alo.ooay)O3 (15). We performed a Rietveld 
refinement (Fig. 1) and provide a crystallo- 
graphic information file (75). 

Next, we describe the hosting diamond and 
then discuss the pressure-temperature (P-T) con- 
ditions of formation of the type davemaoite 
and its composition. The hosting diamond is 
different from diamonds that contain poten- 
tial retrograde products of high-pressure min- 
erals at 0 to 1 GPa remnant pressures (13, 18). 
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Fig. 1. X-ray diffraction pattern and Rietveld refine- 


ment of davemaoite. Davemaoite (dvm), red; (Fe,C), 


olive; wistite (wlis), blue; residual of fit, green. Tick marks 


indicate allowed reflections (same color-coding). The 
weighted refinement factor, Rwp, was 0.046, and the 


structure-factor moduli-based refinement factor, Re, for dvm 


was 0.10. Phase proportions were 35(5), 45(5), and 


15(5)% for dvm, (Fe,C), and wiis, respectively. In dvm, the 
intensity of reflections 100 and 201 indicates sublattice 


disorder with A-cations shifted from site 1b to the 1/8- 











occupied site 8g. This is consistent with accommodation 


of K and Na along with Al and Fe. The paragenesis was 





found in two inclusions. Volumes of analyzed davemaoite 


inclusions ranged from 45.1 to 46.3 A? [equal to 0 to 


6 GPa (7) in pure CaSiO3], but the actual composition of 
davemaoite markedly changes the pressure inferred from 


the volume (9). Coexisting wustite Feo g)Mgo.2a)0 


(22). Indications of weak diffraction around dvm peaks 
near noise level are assigned to minor contributions of 








inconsistent with a tetragonal distortion. a.u., arbitrary 


The unprocessed IR spectrum is shown in fig. S3. Bars 


has a volume of 74.34(1) A°, corresponding to 8 to 9 GPa 
material of slightly larger and smaller volume (probably 
eflecting chemical variation; fit not shown here) but are 


units. (Inset) IR spectrum of davemaoite after subtraction 
of a diamond spectrum collected close to the inclusion. 








The diamond is of type IaAB with frosted 
octahedral faces and trigon features (15). We 
found davemaoite, iron, wiistite, ilmenite, and 
ice-VII in the center of the diamond. Our 
analysis of the N-defect IR bands (fig. S3) (12) 
indicates a low average residence temperature 
(~1500 K) or a short residence time in the man- 
tle, similar to the holo- and cotype diamonds 
of ice-VII (78). Short residence time and low 
average residence temperature are common 
features of lithospheric diamonds, but in sub- 
lithospheric diamonds both parameters act in 
favor of conserving high remnant pressures 
and high-pressure minerals by reducing visco- 
elastic relaxation of the hosting diamond and 
by preventing retrograde transformations. The 
bulk modulus of davemaoite depends on its 
composition (7, 8) and is unknown for the 
given composition. However, coexisting wiistite 
is at a remnant pressure of 8 to 9 GPa (16). 
For a single inclusion of wiistite, this remnant 
pressure would correspond to an entrapment 
pressure of ~40 GPa if pressure were to 
evolve along a purely elastic path (75). How- 
ever, diamond becomes viscoelastic between 
1100 and 1200 K, even at laboratory time scales 
(19). To account for this nonelastic behavior 
of the hosting diamond, we use the method 
of Wang et al. (20), which does not rely on 
initial assumptions about the entrapment 
temperature and uses the P-T paths of sep- 
arate inclusions in the same diamond: wiistite, 
iron, ice-VI, and ilmenite. Using this approach, 
we assessed entrapment conditions of 29 + 
5 GPa at 1400 to 1600 K (15). Because visco- 
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indicate all calculated IR bands. 


elastic processes are path- and time-dependent, 
we cannot exclude a higher entrapment pres- 
sure or temperature. 

We cannot entirely rule out that our chem- 
ical analysis is affected by minor contaminants, 
although we did not observe an XRD signal or 
a marked XRF signal of any phase other than 
davemaoite, wiistite, and iron in the excavated 
region. Furthermore, we note that (i) the low 
Tiis a result unaffected by potential contam- 
ination and (ii) the ?°K signal occurs at the 
same depth as the “*Ca signal of davemaoite. 
We did not observe by XRD alternative 
hosts of K and Ca such as liebermannite and 
harmunite-type (Ca,K,Na)(Al,Si).0, anywhere 
in this diamond. Both of these phases, and any 
phase dominated by K and Ca, give diffraction 
patterns that are very different from the 
perovskite-type pattern we observed. Hence, we 
believe the presence of K and Al in davemaoite is 
not likely the result of a contaminated analysis 
but rather indicates coupled substitution of a 
large and a small cation K,Na + Al,Fe for Ca. 
Generally, a substitution of K for Ca and Al for 
Si shifts the material into the stability field of 
ABOs perovskites with a trend toward high 
crystal symmetry (27). 

Our result indicates that the postspinel phase 
(Ca,Na,K)(AL,Si).O, is not required as a host of 
Ca, alkalis, and Al in at least the upper region 
of the lower mantle. It is possible that type 
davemaoite formed retrograde out of postspinel 
through a reaction (Ca,Na,K)(Al,Fe**,Si).O, + 
Fe° — (Ca,Na,K)(Al,Si)O3 + FeO, but for this 
process one expects the presence of remnant 
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Fig. 2. Reflected light image of the diamond at the beamline. The XRF map 
of CaKa is superimposed. The x and y axes give the sample coordinates in 
micrometers. The red circle indicates the area ablated during LA-ICP-MS 
analysis. The two Ca-rich areas in the circle correspond to diffraction patterns of 
davemaoite. The more intense signal around x,y = 72,-3074 corresponds to 


postspinel in the paragenesis, which is not 
observed. In the deep mantle, davemaoite 
takes on a role similar to that of garnet in the 
upper mantle. Both minerals have a “gar- 
bage can” crystal chemistry that allows them 
to host many elements that are incompatible 
in upper-mantle minerals (6, 7). Our observa- 
tions are fully consistent with the experimen- 
tal results that this mineral dissolves LILEs, 
specifically K (6, 7). Experimental studies that 
were based on peridotite and MORB (mid- 
ocean ridge basalt)-like bulk compositions 
formed davemaoite with lower K content and 
higher Ti content than the type material, which 
is expected for these starting compositions. 
We argue that the low Ti and high K content of 
type davemaoite reflects a different, K-rich 
source composition, possibly resulting from 
deep-mantle metasomatism, which is also 
indicated by the presence of ice-VII and by 
the hosting diamond itself (16). This point 
emphasizes the importance of studying nat- 
ural specimens of high-pressure minerals, 


Tschauner et al., Science 374, 891-894 (2021) 





-100 


because they record a petrologic complexity 
of deep Earth that may not be assessed in ex- 
periments. Depletion of Ti in type davemaoite 
is a possible result of the presence of phases 
that strongly partition Ti, such as liuite, FeTiO;- 
perovskite. Ilmenite has been observed in the 
same diamond at similar remnant pressure as 
the davemaoite-wiistite-iron inclusion, and its 
P-T path intersects the phase boundary of liuite 
(15). Hence, our findings indicate that the 
source rock composition of type davemaoite 
deviated from peridotite and, thus, that chem- 
ical segregation occurs in the lower mantle, 
possibly down to 900 km according to our es- 
timate (fig. $4). This variation in rock composi- 
tion affects heat generation through radioactive 
decay in the lower mantle where davemaoite 
scavenges K, as shown here, and U and Th, as 
experimentally shown (7). 
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Lower mantle “garbage can” 

Calcium silicate perovskite has finally been identified in a natural sample and now has the mineral name davemaoite. 
Tschauner et al. discovered the type mineral trapped at high pressure and temperature as a diamond inclusion (see 
the Perspective by Fei). Structural and chemical analysis of the mineral showed that it is able to host a wide variety 
of elements, not unlike fitting bulky objects into garbage can. Specifically, it has a large amount of trapped potassium. 
Davemaoite can thus host three of the major heat-producing elements (uranium and thorium were previously shown 
experimentally) affecting heat generation in Earth’s lower mantle. —BG 
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Data S1 


Diamond provenance: 


Type locality: Orapa, Botswana. Obtained from David New, originally 81 mg octahedral 
diamond with frosted surfaces and trigon etching features. An outer zone of this diamond 
exhibits a large number of fine growth horizons in octahedral shape, marked by hazy dark zones 
of few um thickness. The kernel of the sample is homogeneous both optically and in the XRF 
map, besides the inclusions shown in Fig 2 and Fig S1. The examined sample was a 0.5 mm 


thick platelet that had been laser-cut and polished along the (001) base plane. 
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Fig S1, Top: XRF Fe- (left) and Ti- (right) Ka. Both inclusions of davemaoite contain Fe 
(mostly as iron) and very little Ti. Iron is the most common inclusion in this diamond, consistent 
with its formation in the regime of Fe** autodissocation and redox freezing (23). Fe-signal at 
lower right corner is surficial contamination (visible in the microscope image in Fig 2) that was 
blasted away upon pre-ablation cleaning prior to LA-ICP-MS analysis. Based on XRD the 
elongated feature in the upper left corner around Y = -3000 and -170 < X<-120 is from ilmenite 


and the feature centered around Y ~ -3030 and X ~ -140 is from taenite. 


Bottom: Time-resolved signals of *°Fe (red) and *?Cr (black) (Left), and *°K (black) and “Ca 


(red) (Right). K and Ca signal occur at the same depth. 


Infrared Spectroscopy: 


Mode analysis of davemaoite gives three triple-degenerate IR-active modes T1u, one tripe- 


degenerate silent mode T2u, and no Raman-active modes based on the SAM algorithm (24): Topt = 


Tou + 3T iu = all twelve optical modes (the five atoms in the unit cell give 15 modes, whereof three 


modes are acoustic). 


We conducted IR-absorption at the locations where davemaoite occurs and the 
surrounding matrix for comparison (Fig S2). Infrared (IR) spectra were collected in transmission 
at beamline 1.4 at the Advanced Light Source, Lawrence Berkeley National Laboratory, with a 
Nicolet Magna 760 FTIR bench and a Nic-Plan IR microscope with a 32x magnification 
Schwarzschild objective, with 1 cm‘! resolution and a HgCdTe detector with a KBr beamsplitter. 
Apertures were set to 20x20 um? spatial resolution in order to maximize spatial resolution but 
obtain signal in the spectral range above 500 cm’. A synthetic type II diamond was used for 
assessing the absorbance. Consistently with mode analysis, no Raman peaks were observed at the 
location of davemaoite (Fig S2). Raman-spectra were collected with a Renishaw M1000 Micro 
Raman Spectrometer System at 514 and 783 nm at Caltech and with a custom-built confocal micro- 


Raman system at GSECARS with 473 nm excitation wavelength. 


The most recent ab initio calculations of phonons of the synthetic analogue of davemaoite, 
CaSiO3-perovskite, were conducted by Caracas et al. (/8) and Sun et al. (25), the latter for high 
temperature. In both studies the phonons of the low-symmetric space-group variants of CaSiO3- 
perovskite were calculated and no mode assignments were given. Therefore, we reproduced all 
reported mode energies with a force-constant model, applied it to the cubic structure and identified 
the phonons at 270, 740, and 890 cm”! at I that were calculated by Caracas et al. (18) with the 
three T:u IR-active modes of cubic CaSiO3-perovskite. These modes correspond to Ca-sublattice 
breathing, bending of O- relative to Ca-sublattice, and to rotational bending of the Si-O- bonds 
with simultaneous distortion of the Ca-sublattice, respectively. The absorption band at 740 cm’! is 


expected to be most intense based on the phonon density of state (25). Within the accessible spectra 


range we observe the band at 740 cm’! and see an indication for presence of the band at 890 cm’! 
(Fig 1). The mode at 270 cm’! cannot be measured because the necessary diffraction limit would 
require degrading the spatial resolution to ~ 100 um and the inclusions are smaller than the 


diffraction limit at that wavelength. 


Nitrogen defects: Absorbance between 1000 and 1500 cm’! was assessed through 
measurements of a synthetic type II diamond. Spectra were scaled to an equal thickness of 10 mm 
after background subtraction. The contributions of A-, B-, and D-type defect bands were assessed 
based on the defect type spectra from Taylor et al. (26). The amount of nitrogen in A- and B-type 
defects were estimated based on the calibrations by Boyd et al. (27, 28). The amount of N in D- 
type defects (29) was within uncertainties of the amounts of A- and B-defects of 93411 and 52+6 
ppm, respectively. An average mantle residence temperature of the diamond was estimated based 
on the nitrogen contents from A- and B-type defects and the corresponding percentage of A-defects 


and assuming an age of the diamond of between 0.8 and 1.6 Ga using the rate equation in (26). 


Analysis of chemical composition 


Chemical compositions of the diamond inclusions were measured using a Thermo Element 
XR™ ICP-MS coupled with an Elemental Scientific Lasers (ESL) New Wave™ UP193FX 
excimer (193 nm) laser ablation system at the Plasma Analytical Facility of the National High 
Magnetic Field Laboratory, Florida State University. The laser ablation system was operated using 
spot mode, with a spot size of 100 um, 50 Hz repetition rate, and 80% energy output corresponding 


to 12.6 Jecm”. All peaks were monitored in medium mass resolution (m/Am= 4,000), so that the 


complicated carbon- and argon-related molecular interferences on low mass isotopes were 
resolved. Jet-sampler and X-skimmer cones were used to compensate for the diminished signals 
in medium mass resolution. Peaks of Li, "Be, '7C, °C, '°O, “Na, 4Mg, 7’Al, *8Si, 3'P, °°S, °°K, 
4Ca “Se, 7Ti, 5'V, 2Cr, Mn, Fe, °Co, Ni, Cu, Zn, Ga, “Ge, “Rb, Sr, °Y, Zr, Nb, 
I3Cg B8Ba, La 49Ce, 4Ipp, 45Nd, 47gm, !3 Eu, Tb, !Gd, !Dy, Ho, Er, 'Tm, '#Yb, 
Tu, 8°HF, '8!Ta, 2°8Pb, 7?°Th, and 8U were monitored using magnet peak jumping technique 
with the Element’s EScan mode by varying the acceleration voltage. Peaks were monitored using 
10% mass window with an acquisition time of 0.1 second per peak. *°Ar*°Ar peak was used as the 
lock mass. Blanks were measured by ablating the same diamond in an inclusion-free spot, as well 
as on a piece of broken diamond anvil. The blank corrections on most major elements are on the 
order of percent level or lower. The measured count ratios referenced to *°Si were converted to 
elemental abundance ratios referenced to Si using glass standard GSE-1G (30). The elemental 
compositions were then calculated assuming that the sum of all the major oxides is 100% (3/). 
During our LA-ICP-MS measurements, !°O*/4?Ar*°Ar* was about 3 and ?C*/°Ar*?Ar* was about 
0.7 when diamond was ablated. The high '°O*/*’Ar*?Ar" ratio is due to air entrainment. *°Ar* was 
not monitored, but *’Ar*’Ar*/*°Ar" is typically less than 0.001. Consequently, during our LA-ICP- 
MS measurements, Ar* was the dominant species in the plasma whether silicate samples (O-rich 
matrix) or diamond samples (C-rich matrix) were ablated. All elements measured in our study have 
first ionization potentials (FIP) < 10 eV. Using the Saha Equation, we found that they are nearly 
completely (>90%) ionized in the Ar plasma regardless of their matrix, silicate or carbon. 


Therefore, standardization based on a silicate standard does not cause systematic errors. 


The effect of adding carbon to an ICP Ar plasma was studied by Fliegel et al. (32), who 


reported an overall increased sensitivity and, specifically, that the sensitivities of low FIP (<10 eV) 


elements uniformly enhanced by a factor of several, and those of high FIP (+10 eV) by up to more 
than one order of magnitude. We noted that Fliegel et al. (32) added 1.0 + 0.4 ml/min methane in 
their experiments, which is equivalent to a carbon flux of 5 + 2 x 10° ng/min. For comparison, 
based on the size of the ablation pit on the diamond, we estimated a carbon flux of several hundreds 
of ng/min during our LA-ICP-MS measurements. If we scale the sensitivity enhancement observed 
in Fliegel et al. (32) to our measurement, we expect to see at most several percent enhancement in 
the high-FIP elements, which is within our analytical uncertainty of several percent. In addition, 
most elements reported in this study (Table S1) have low FIP. Addition of carbon might uniformly 
enhance their sensitivity. However, this effect, if any, is canceled out during our data reduction 


process, because we normalize all analyzed elements to Si. 


In summary, it is justified to use a silicate standard, GSE-1G, to calibrate the compositions 


in diamond inclusions, since the effects of adding the small quantity of carbon are negligible. 


The composition of the individual minerals was obtained from the composition of the bulk 
inclusion by a) the phase proportions assessed by Rietveld refinement (33), b) varying the Fe- and 
Mg-content of wiistite within the range given by the Rietveld refinement of (Feo.s1),Mgo.271))O, c) 
the analytical uncertainties, and d) charge balance of davemaoite. The at-% in Table S1 are given 
as percent of the inclusion total composition. The observed tetragonal distortion of the unit cell of 
iron indicates that it contains a few at% C (34). The uncertainties in phase proportion are the largest 
source of uncertainty, but because of the additional constraint of charge balance of davemaoite, 


they mostly affect the composition of iron and wiistite. 


Table S1: Composition of davemaoite, iron and wiistite. Composition is given in at% of 


the total composition of the two inclusions in average. 



































at-% Davemaoite | Iron Wiistite 
Ca 18.27(53) - - 

K 8.71(30) - - 
Na 2.6(2) - - 

Fe 4.7(1) 9.5(9) 1.2(3) 
Al 3.32(17) : ? 
Mg 2.43(7) - 0.15(5) 
Cr 1.8(8) 0.2(1) - 

Ti 0.30(3) - ? 
Mn 0.23(3) - > 0.03 
Si 42.2(8) 0.3(2) - 

















The empirical formula on the basis of 3 O pfu is (Cao.43(1)Ko.2001)Nao.06F e0.11(1)Alo.osMgo.o6Cro.04(2)) 


(Si1.0(2)Alo.001))O3 . There is no ABO3 perovskite with the observed unit cell dimensions except 
for (Ca,...)Si03. To our best knowledge GeCFe3, GeCCos, and GaCNi, are the only perovskites 
with similar cell dimensions, but Ga and Ge in the present sample were below their detection 


limits and none of these synthetic phases contain Ca. 


Structure analysis: 





Single-crystal X-ray studies could not be carried out owing to the polycrystalline nature of the 
holotype material and its sub-um-scale grain size. X-ray powder diffraction data of the three 
inclusions were collected with synchrotron radiation at beamline 34-IDE at the Advanced Photon 
Source, Argonne National Laboratory with focused primary beams of 0.5166 A wavelength and 
a MARI165 CCD area detector, calibrated, corrected for geometric distortion with Dioptas (35) 
and subsequently integrated with Dioptas. The specimen was mounted on a motorized x-y-z 
stage with sub-tum motion precision and centered in the beamwaist of the focused X-ray beam 
based on a high-magnification optical system that was aligned previously with an X-ray 
fluorescent YAG crystal. The X-ray beam waist is ~ 300 um in depth and samples almost the 
entire thickness of the sample without significant loss of X-ray flux density. Diffraction data 
were collected in mapping mode over a 40x40 tm? grid in 2 um steps with 10s acquisition time 
per frame. Frames from adjacent diamond matrix were used for diffraction frame background 
subtraction. The three inclusions gave diffraction patterns that matched an XRD pattern of the 


cubic perovskite-type. 


The structure of synthetic perovskite-type CaSiO3 is well established (9) and the phase was 
easily identified. In the perovskite-type CaSiO3 structure, all atoms reside on special sites (Table 
S2). No indication for spontaneous structural distortion was observed, consistent with absence of 
observable Raman-shifts. Cubic perovskite-type CaSiO3 has no Raman-active modes. Refined 


parameters were background, pseudo-Voigt profiles, Biso, and preferred orientation along 111. 


Rietveld refinement converged to Rwp = 0.046 and Rp = 0.035 with a x” of 17.05 for 2041 
observations. The final Rr was 0.098. Profile parameters were 1.10(2), 0.05(1), 0.014(1) for the 
Cagliotti terms, and the Lorenzian mixing term was 0.485(5). These values are added in the 


comment lines of the cif file. We obtained the following phase proportions (in vol-%): 


davemaoite 30-40%, iron (Fe,C,Si,Cr) 40-50%, wiistite 10-20%. Preferred orientation along 111 
was fitted with the March-Dollase approach (orientation factor 0.60(5)). A Le Bail extraction 
converged to Rp = 0.0378 and the extracted |F(hkl)| were used for assessing Rr to 0.098 (Table 
S2). LeBail extraction and Rietveld refinement were conducted with PowderCell (33) and GSAS 
(36). Figure 2 shows presence of additional weak diffraction signal around davemaoite- peaks 
and we examined the possibility of structural distortion of this phase. Two distinct paths were 
chosen: 1) We used the extracted |F(hk1)| for reversed Monte Carlo modeling without crystal 
symmetry bias and without potentials. We conducted both, local and global optimization (37) for 
sampling and extended configurational space and in both cases we recover the cubic structure. 11) 
We forward modelled tetragonal, orthorhombic, and rhombohedral cells. The former two were 
modeled based on the work by Chen et al (9) and Kurashina et al. (38). For the latter we mapped 
the cubic cell onto a rhombohedral subgroup. Then we conducted Rietveld refinements for each 
of these model structures. However, in each of these refinements the fractional atomic 
coordinates converged toward the cubic arrangement or, if constrained, predicted non-observed 
peak intensities that could not be modeled by preferred orientation. The wRy did not improve for 
these model structures. The weak additional signal is best modeled by the presence of an amount 
of < Svol% of cubic davemaoite of ~ 2.5% larger or smaller cell volume which may reflect minor 
variations in composition or local stresses for small, separate inclusions in the vicinity of the 
larger inclusions. We did not include the refinement of this additional minor phase in the 
Rietveld fit that is shown in Figure 2 and also not in the extraction of |F(hkl)| in table S2 and the 
cif-file. WRp and Rr slightly improve with this addition but the improvement is statistically 


marginal while it adds further refineable parameters (profiles, cell volume, and phase 


proportions). Therefore, we did not include this in the final results here. Observed and calculated 


normalized structure factor moduli of extracted reflections and their d-spacings are listed in 


Table S2 (a) Miller indices, d-spacings, observed and calculated structure factor moduli of 


davemaoite type-material (a), iron (b), and wiistite (c). M= multiplicity, calculated |F(hkl)| are 


normalized relative to the observed |F(200)| = 1000. The final RF = 0.098. The d-spacings 


represent the centroid positions of the peaks in the LeBail fit. The uncertainties of the d-spacings 


are < 0.001A. 


d(A) |Fcalc| 


3.583 


2.534 


2.070 


1.791 


1.603 


1.463 


1.267 


1.194 


1.194 


|Fobs|_ 1-sigma(Fobs) 
Dio) Qe “Bs 

631.6 537.85 93.75 

635.75 726.65 -90.9 

1152.65 1000 152.65 

19.2 18.45 0.75 
463.55 426.4 37.15 
770.75 755.6 15.15 
15.4 


15.35 0.05 


15.4 15.35 0.05 


Table S2 (b) List of observed and calculated |F(hk])| for coexisting wiistite (Fe0.8(1)Mg0.2(1))O, 


partial site occupancy of Fe was refined. The unit cell parameter a = 4.202(1)A. Rr = 0.091. 
h k l M d(A) |Fcalc| — |Fobs|_ 1-sigma(Fobs) 

1 1 1 8 2.426 564.3 475.6 88.7 

2 0 0 6 2.101 1066.8 1000.0 66.8 

2 2 0 12 1.485 792.1 811.5 -19.4 

3 1 1 24 1.267 372.7 429.0 -56.3 


2 2 ye) 8 1.213 625.2 705.0 -79.9 


Table S2(c): List of observed and calculated |F(hkl)| for coexisting (Fe,C,Si,Cr), iron-O, Im m 


m, RF = 0.098. The unit cell parameters are a = 2.904(4), b= 2.892(4), c = 2.863(4) A. 
1 1 0 4 2.049 993.7 933.9 59.8 

1 0 1 4 2.039 855.6 851.9 3.7 

0 1 1 4 2.035 1000.1 1000.0 0.1 

2 0 0 2 1.452 553.2 481.5 71.7 

0 2 0 2 1.446 848.7 798.9 49.8 

0 0 2 2 1.432 840.1 888.9 -48.8 


2 1 1 8 1.182 369.0 381.0 -12.0 


1 2 1 8 1.179 528.1 579.4 -51.3 


1 1 2 8 1.173 585.6 658.7 -73.1 
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Fig $2: Raman spectrum taken at the location of davemaoite. Only the Raman shift of diamond 


is visible. 
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Figure S3: Infrared spectrum of N-defect related absorption bands. Blue: A-type defects, red: B- 
type defects, green: D-type defects. Right panel: Spectrum of diamond at location of davemaoite 
(black) and reference spectrum taken 50 micrometer afar (red). The spectra were offset and 

rescaled. The arrow indicates the principal absorption band of davemaoite that is not observed in 


the reference spectrum. 


Assessment of the entrapment pressure and -temperature of the diamond inclusions 
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Fig S4: Estimate of entrapment conditions based on Wang et al. (2/). Blue: ice-VII and fluid 
H20, olive: wiistite, red: ilmenite and liuite (calculated, /8), purple: bec- and fec-iron. Up to the 
elastic-viscoelastic transition in diamond (grey line, linear interpolation of yield strength at 10 
GPa confinement pressure, 20) the paths are isomekes (37) with foot pressure at 300 K given by 
the measured volumes and experimental isotherms (2/) after correction for elastic relaxation of 
the hosting diamond (39). Once viscoelastic yielding of diamond becomes noticeable at the stress 
of the inclusions their actual P-T paths are time(rate)-dependent. We used adiabatic paths as 
lower limits representing complete elastic equilibration of inclusion, hosting diamond, and 
surrounding rock (2/). The adiabats intersect the average residence temperature of the hosting 
diamond (hachured area) of 1500 + 100 K at 25 GPa, 1600 K and 35 GPa, 1400 K which values 


define a lower limit of the entrapment conditions where the actual paths of the different 


inclusions converged (red ellipse in Figure S4). However, it cannot be excluded that the actual 
paths could have been steeper than adiabatic or entrapment occurred at higher than average 


residence temperature. 
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